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Abstract 

The p_L-distributions for the dimuon production at the LHC are 
calculated in the Randall-Sundrum scenario with a small curvature k. 
The pp collisions at 7 TeV and 14 TeV are considered. The widths 
of massive graviton excitations are taken into account. It is found 
that the LHC discovery limit on 5-dimensional gravity scale M5 is 
equal to 5.5 TeV for -y/s = 7 TeV and integrated luminositiy 5 fb~^. 
For = 14 TeV and integrated luminosities 30 fb^^ and 100 fb^^, 
the search limits in the dimuon events are 12.0 TeV and 14.6 TeV, 
respectively. Contrary to the standard RS model, these limits are 
independent of provided k <C M5. 

1 Extra dimension with the small curvature 

The goal of the present paper is to estimate gravity effects in the dimuon 
production at the LHC in a scheme with a warp extra dimension, and to 
derive search limits on a fundamental 5-dimensional gravity scale M5. 

The standard Randall-Sundrum (RS) model [1] is the theory with one 
extra dimension (ED) in a slice of the AdSs space-time. It predicts a series 
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of massive Kaluza-Klein (KK) resonances with the hghtest one, mi, around 
one TeV. A lot of efforts were made in order to find effects coming from 
warped ED. The bounds on 5-dimensional Planck scale M5 and/or mi were 
obtained both at the Tevatron ^ and recently at the LHC In particular, 
high-mass dilepton [3J, diphoton [3] and tt resonance events [5J were searched 
for in pp collisions at ^/s = 7 TeV with the ATLAS and CMS experiments 
(see also [H], [7]). 

In the present paper we will study the RS scenario with the 5-dimensional 
Planck scale M5 in the TeV region and small curvature [8]-|12]l^ 

K < M5 . (1) 

In such a scheme, the background warped metric looks like 

ds^ = e^<^'^-\v\) dx" dx" + dy^ . (2) 

Here y = r^O (— tt < < tt) is the 5-th dimension coordinate, rc is the size 
of the ED, Tj^y is the Minkowski metric. The points {x^,y) and {x^,—y) 
are identified, so one gets the orbifold S^/Z2. The parameter k, defines the 
5-dimensional scalar curvature of the AdSs space. In what follows, we will 
call it curvature. 

We are interested in the RS scheme with two 3D branes with equal and 
opposite tensions located in 5-th dimension at the points y = nrc (called TeV 
brane) and y = (called Plank brane). The SM fields are constrained to the 
TeV brane, while the gravity propagates in all five spatial dimensions. 

It is necessary to note that metric is chosen in such a way that 4- 
dimensional coordinates x^ are Galilean on the TeV brane where all the SM 
fields live, since the warp factor is equal to unity at y = vrrc. Due to the warp 
factor in the metric (|2]), the gravity is strong on the Planck brane, while it 
is weak on the TeV brane. 

By integrating 5-dimensional action in variable y, one can get an effective 
4-dimensional action which, in its turn, leads to so-called hierarchy relation 
between the 5-dimensional reduced gravity scale M5 and reduced Planck mass 
Mpi: 

Ml, = (e'"""^ - 1) . (3) 



^The ratio of the S-dimensional scalar curvatures in our scheme TZ to the scalar curva- 
ture in the standard RS scheme TZq is given by TZ/TZo ~ {k/M^Y ^ 1- 
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The fundamental gravity scale M5 (i.e., Planck scale in five dimensions) and 
reduced scale M5 are related as follows: 

M5 = (27r)^/2M5 ~ 1.84M5 . (4) 

In order the hierarchy relation (|3]) to be satisfied, it is enough to take r^K ~ 
8-^9.5, that corresponds to ~ 0.15-i-1.8 fm. In what follows, we choose M5 
to be of the order of several TeV to tens TeV, while k is allowed to vary from 
hundred MeV to tens GeV. Thus, no new parameters of the order of Mpi (as 
in the standard RS model [1], in which k ~ M5 ~ Mpi) are introduced in our 
scheme. 

From the point of view of a 4-dimensional observer located on the TeV 
brane, in addition to the massless graviton, there exists an infinite number 
of its Kaluza-Klein (KK) excitations, gIi""*, with the masses 

m„ = n = 1,2. . . , (5) 

where Xn are zeros of the Bessel function Ji{x), and n is the KK- number. 

The interaction of the KK gravitons with the the SM fields on the TeV 
brane is described by the Lagrangian 

1 °° 
Mpi ^ K 

were T'^^ is the energy-momentum tensor of the SM fields. The parameter 

can be regarded as a physical scale on the TeV brane where all SM fields live. 

The scalar field, radion, is omitted in the Lagrangian ([6]). In our scheme, 
it practically decouples from the SM fields since its coupling on the TeV 
brane is proportional to [8j 

In particular, for k = 1 GeV and M5 = 1 TeV, one gets (^rad ~ 1/(100 TeV). 

As for processes with a real or virtual production of the KK gravitons, the 
smallness of their coupling to the SM particles, l/Ajr, is compensated by a 
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large number of real gravitons or by infinite tower of virtual KK excitations. 
As a result, a magnitude of the cross sections is defined by the scale M5, not 
the scale [8] (as an illustration, see fl2T]) ). 

As it was already mentioned above, the 4-dimensional coordinates are 
Galilean on the TeV brane. Thus, a correct determination of the masses 
on this brane can be achieved [13]. In particular, observed masses coincide 
with their Lagrangian values and do not depend on a coordinate rescaling, 
if covariant equations and invariant distances are used |14] . 

The RS-like model under consideration has other interesting features. For 
instance, the spectrum of the KK gravitons ([5]) is very similar to that in the 
ADD model with one ED [T3]. Moreover, all matrix elements can be formally 
obtained from corresponding matrix elements calculated in the ADD model 
with one flat extra dimension by using following replacements [H [TO] : 

M4+1 ^ (27r)-i/3 M5 , R^^ (ttk)-' . (9) 

Here M4+1 is a 5-dimensional reduced Planck scale, Rc being the radius of 
the compact fiat dimension. 



2 Gravitons in dimuon production at the LHC 

In this section we will study aproduction of two muons with high transverse 
momenta in pp-collisions [16] □ 

pp^H~^H~ + X. (10) 

We are interested in a possible excess in p^-distribution of the final leptons 
coming from gravity interactions in 5-dimensional space-time with the small 
curvature. The theoretical framework was described in details in the previous 
section. 

The dilepton production is a favourable process to look for gravity effects 
predicted by theories based on a space-time with more than four dimensions. 
The Tevatron data on p^-distribution in dilepton production already pro- 
vided us with bounds on parameters of these theories [17j. As for dilepton 
production at the LHC, the search for EDs resulted in new limits on so-called 
large EDs [18]. However, most efforts of ATLAS and CMS Collaborations 

■^The production of a lepton pair in pp (or pp) collision is often called Drell-Yan process. 
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were made to search for massive resonances in such events predicted by the- 
ories with the warped metric (see, for instance, [SIE])- The phenomenology 
of the RS model with the resonant KK spectrum in the TeV region can be 
found in [19]. The determination of the spin-2 RS gravitons against spin-1 
and spin-0 resonances in lepton-pair production was done in |20] . 

The goal of this paper is to estimate gravity effects in the dimuon pro- 
duction at the LHC in the RS scenario with the small curvature [H]-[in]- The 
differential cross section of this process ( ITO|) is given by 



da 
dp± 



{pp ^ /.V- +X)=2p^Yl I / — 



Xi] 



X /fe/p(yU^^/a;l) ^{ab ) , (11) 

at 

with the transverse energy of the muon pair equals to 2pj_. Here fa/pifJ''^, xi) 
is the distribution of the parton of the type a in momentum fraction Xi 
inside the proton taken at the scale fi. da/dt denotes the cross section of 
the partonic subprocess ab — )■ which is described by the Mandelstam 

variables s, t and u (s + t + m = 0). 

We introduced two dimensionless quantities, 

x±_ 



T = XiX2, (12) 

where X2 is the momentum fraction of the parton b in (fTTj) . Without cuts, 
integration variables in ( ITTl) vary within the following limits 

x'i<T <1 , 

T <Xi<l . (13) 

After imposing kinematical cut, the integration region becomes more com- 
plicated (see Appendix A). 

The contribution of the virtual gravitons to lepton pair production (/ = 
e or /i) comes from the quark- ant iquark annihilation, qq G*-"-' — 
and gluon-gluon fusion, gg ^ G*^"-* — )■ l^l". The corresponding partonic 
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cross sections are (see, for instance, Appendix A in [5]) 



Z^^^"^^ ^"^ = ^^^^^i^ ' 

da ^ , t(s + t)(s^ + 2st + 2P) , ^,^,,2 , 

(gg^tr) = ^-S{s)\ 14 

dt Zbbns^ 

where 

1 ~ 1 

'^^"^ = Af g 7^<T7^ ^^^^ 

is the invariant part of the partonic matrix elements, with r„ being total 
width of the graviton with the KK number n and mass m„ |TU] : 

in 

\ = Vmn{^] , V-Om. (16) 

Note, because of the universal coupling of the KK gravitons to all SM fields, 
the function S{s) is the same for all processes mediated by s-channel virtual 
gravitons. 

In the RS scenario with the small curvature, sum ([T5l) was calculated 
analytically in [10]. At s ~ M5 ^ k, it looks like 

1 sin 2yl + 2 sinh 2e 
^ ~4M|7i cosM + sinh^e ' ^ ^ 

where 

A=^, e = ^-(^y. (18) 

It is interesting to note that the magnitude of S{s) is defined by the 
fundamental gravity scale M5, while the scale A^, which is contained in the 
Lagrangian ([6]), did not appear in f lT7|) . Let us underline that only for e 1 
we come to the expression which can be obtained by neglecting widths of the 
massive gravitons [zero width approximation)^ 

liaSis)- J-^ , Re5(s)~0. (19) 



^In fact, zero width approximation is justified if inequality s > 3.5 M5 is satisfied |10) . 
Under this condition, Re5/Im5 < 0.05. 
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In what follows, we will consider M5 > 2(4) TeV for the LHC energy ^/s = 
7(14) TeV. In such the parameter e in (ITSl) obeys inequality e < 2, 

and one must use the exact formulae f|T7|) . f|T8|) instead of the approximate 
expressions f|T9l) FI 

As one can see from f|T7|) . the graviton spectrum in our model consists 
of very sharp resonances whose widths are proportional to (/t/Ms)^. The 
contribution from one resonance to the p_L-distribution can be estimated as 
follows (more details can be found in \T2\): 



(icT (grav) 



dpA 



(20) 



Since the total number of graviton resonances which contribute to the cross 
section is proportional to y/s/n, we obtain at fixed x± 

da{gmv) 1 

^ Ml ■ 

Let us stress that in our scheme the gravity cross sections do not depend 
on the curvature k, provided k <^ Msj^l in contrast to the RS model with the 
large curvature in which n ~ M5 ~ Mpi p.j. 

In order to obtain search limits for the LHC, we calculate contributions 
from s-channel graviton resonances to the p^-distributions of the final muons 
for different values of 5-dimensional Planck scale M5. We use the MSTW 2008 
NNLO parton distributions [21], and convolute them with the partonic cross 
sections f lT^ in flTT]) . The PDF scale is taken to be equal to the invariant 
mass of the muon pair, /i = Mi+i- = a/J. 

We also impose the cut on the lepton pseudorapidities which is used by 
ATLAS and CMS Collaborations in selecting dimuon events [6]: 

1^1 < 2.4. (22) 

The limits on variations of the variables r and Xi in (fTT]) . resulting from this 
cut, are derived in Appendix A. 

The reconstruction efficiency of 85% is assumed for the muon events 
Note that this assumption is not very crucial for our scheme. For instance. 



^ Since a mean value of the partonic energy \/l is less than the collision energy ^/s, an 
effective value of the parameter e is even less than 1. 
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Up to very small corrections of the type 0(k/M; 
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given the efficiency varies from 85% to 80%, the LHC search hmits on M5 
diminish by 1% only, due to the strong power-hke dependence of the gravity 
cross section on M5 fl2T]) l^ 

In Fig. [1] and Fig. [2] we present the resuhs of our calculations of gravity 
cross sections for the dimuon production at the LHC. The gravity mediated 
contributions to the cross sections do not include the SM contribution. The 
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Figure 1: The KK graviton contribution to the dimuon production at the 
LHC for several values of 5-dimensional reduced Planck scale (solid curves) 
vs. SM contribution (dashed curve) at i/i = 7 TeV. 



ratios of the gravity induced cross sections to the SM one is shown in Fig. |3l 

The differential cross section of the process under consideration is repre- 
sented in the form: 

da = da{SM) + rf(T(grav) + rfcT(SM-grav) , (23) 

where the last term comes from the interference between the SM and graviton 
interactions. Since the SM amplitude is pure real, while the real part of each 
graviton resonance is antisymmetric with respect to its central point, the 
interference term in (1231) has appeared to be negligible in comparison with 

^Namely, a significant variation of da{gTa,v)/dp± can be achieved by a relatively small 
variation of M5. 
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Figure 2: The same as in Fig. [H but for the energy ^/s = 14 TeV. 
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Figure 3: The ratio of the gravity induced cross sections to the SM dimuon 
cross section at the LHC as a function of the muon transverse momentum 
for several values of 5- dimensional reduced Planck scale. 
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the pure gravity contribution (the second term in 0231) ) after integration in 
partonic momenta. More details can be found in [12j. 

The next Fig. H] demonstrates us that an ignorance of the graviton widths 
would be a rough approximation. As one can see, it results in very large 
suppression of the cross sections. The reason of this suppression lies partially 
in the fact that 



(i(T(grav) 
dp± 

while in zero width approximation 

(icr(grav) 



1 

pI 



Ms 



dp± 



zero width 5 



M5 



as one can derive it from formulae (fTTj) . and (fT9l) . 
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Figure 4: The contributions to the dimuon events from the KK gravitons with 
nonzero widths (solid curves) vs. contributions from zero width gravitons 
(multiplied by 10'^, dashed curves) for the different values of the reduced 
5-dimensional Planck scale. 



To take into account higher order contributions, we use a A'-factor 1.5 
for the SM background|3 while a conservative value of -ft" = 1 is taken for 

^For the dimuon channel, the expected SM background consists of Zf-f* , tt and diboson 
events {Z/j* being dominant), while {W + jets) and QCD background is negligible [23) . 



10 



the signal. Let Ns{Nb) be a number of signal (background) dimuon events 
with p± > p]^*. Then we define the statistical significance S = Ns / a/ Nb , 
and require a 5a effect. In Fig. |5] the statistical significance is shown for 
■\/s = 7 TeV as a function of the transverse momentum cut p™* and reduced 
5-dimensional gravity scale M5. Figs. [6] and [7] represent S for the dimuon 
events at y/s = 14 TeV. 




, ItV 



Figure 5: The statistical significance S for the dimuon production at the 
LHC for y/s = 7 TeV and integrated luminosity 5 fb~^ as a function of the 
transverse momentum cut p™^ and reduced 5-dimensional gravity scale M5. 
The plane S" = 5 is also shown. 



The aurthors of Ref. [9] proposed a search in the dilepton channel by 
cutting on the invariant masses Mi+i- above 2 TeV. They found a 5 sigma 
bound on M5 = 4 TeV (without taking into account finite graviton widths). 
In our case, we have the lower bound on Mi+i- = 2p± and integrate in 
variable r = Mf+^^/s. 

3 Conclusions and discussions 

In the present paper the RS-like scenario with the small curvature of the 
space-time [8]- [10] is considered. In such a scheme, the reduced 5-dimensional 
Planck scale M5 can vary from few TeV to tens TeV, while the curvature k 
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400 7 



Figure 6: The same as in Fig. [5l but for ^/s = 14 TeV and integrated 
luminosity 30 fb~^. 




Figure 7: The same as in Fig. [5l but for a/s = 14 TeV and integrated 
luminosity 100 fb~^. 
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is allowed to vary from hundred MeV to few GeV (in fact, the only condition 
K <^ M5 should be satisfied). The mass spectrum and experimental signature 
of the model are similar to those in the ADD model [15] with one flat extra 
dimension. 

The p^-distributions for the muon pairs production with high p± at the 
LHC are calculated for the collision energies ^/s = 7 TeV and y/s = 14 TeV 
(see Figs. [HIS]). The importance of the account of the KK graviton widths 
is demostrated (see Fig. Hj). The statistical significance as a function of 
the reduced 5-dimensional Planck scale M5 and cut on the lepton transverse 
momentum p™* is calculated (see Figs. [MZl)- 

By using these results, we can obtain the following discovery limit for the 
7 TeV LHC and integrated luminosity of 5 fb~^ in the dimuon production: 



In deriving ([26]) and ([27]), we used the relation M5 = {27rY/^M5 

It is important that these bounds on M5 do not depend on the curvature 
K, contrary to the standard RS model [1] in which estimated bounds on M5 
depend significantly on the value of the ratio k/Mpi. 

Previously, analogous bounds were obtained for the diphoton production 
at the Tevatron and 14 TeV LHC [12j. Recently, results on the dilepton 
production at very high luminosity (HL-LHC) were presented [21]. 

Our approach can be directly applied to the production of the electron 
pairs at the LHC. The only thing to do is to take into account that for electron 
samples the transition region 1.37 < \ri\ < 1.52 (1.44 < \ri\ < 1.57) between 
the ECAL barrel and endcap calorimeters is usually excluded in the ATLAS 
(CMS) experiment, while the cuts |?7| < 2.47 and \ri\ < 2.5 are imposed by 
the ATLAS and CMS experiments, respectively (see, for instance, [6]). 
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M5 = 5.5 TeV . 



(26) 



Correspondingly, we obtain for the 14 TeV LHC: 




(27) 
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Appendix A 



In this section we will calculate the integration region in ( [TT|) in a case when 
a cut on final lepton (pseudo) rapidities is imposed. Let us consider a pro- 
duction of a lepton pair with high transverse momenta p± in pp collision: 



pp ^ I + X , / = e, /i . 
In the c.m.s. of the initial protons, their 4-momenta are 

Pi2 = ^(1, 0, ±1) , 



(A.l) 



(A.2) 



where y/s is the collision energy. Let us call this system lab-system, contrary 
to the c.m.s. of the final leptons. In what follows, we neglect the masses 
of the protons and leptons with respect to ^/s, as well as inner transverse 
momenta of the partons inside the protons with respect to p±. 

Let Xi, X2 be momentum fractions of the partons gi, q2 inside the colliding 
protons. Then the 4-momenta of the final leptons are: 



'1,2 



(A.3) 



with 



,11 

,0 _ 
*i,2 - 4 



(Xi - X2) ± (Xi + X2)\l 1 - 



(Xi + X2) ± {Xi - X2)\/ 1 — 



X I 



(A.4) 



The dimensional variables r and x±_ were defined in the main text f|T2l) . They 
obey kinematical inequalities f lT3|) . The invariant mass of the lepton pair is 

S = S XiX2- 

In the c.m.s. of the final leptons (herein called cm- system), the proton 
momenta look like 



Pi 



P2 



2 



Xi 
Xi 
X2 




Xi 



, 0, 



IXi 
X2 



(A.5) 



16 



while the parton momenta are 



qi,2 = {^/XlX2, 0, ± yJXiX2) • 

Correspondingly, the momenta of the leptons looks like 



'■1,2 



ST 



XI 
T 



(A.6) 



(A.7) 



Note, since the lab-system and cm-system are related by a Lorentz boost 
along the beam axis, the partons momentum fractions remains the same in 
both systems. 

Let us turn to the lab-system. The rapidity of the lepton pair in this 
system is equal to (1/2) ln(a;i/a;2). This means that the lepton rapidities are 
given by 



yi 



2/2 




(A.8) 



We are interested in the production of the leptons with the high momenta 
(p± > 200 GeV). That is why, in order to estimate the integration region in 
f lTTj) . one can safely impose a cut on lepton rapidities yi^2, instead of using 
cut on their pseudorapidities^ 



I 2/1,2 |< Vcut ■ 

This bound is equivalent to the following inequalities: 



T 



where 



A 



x± exp{y, 



cut J 



(A.9) 



(A.IO) 



(A.ll) 



^Remember that the rapidity is defined as y = (1/2) \n[{E + pz)/{E — Pz)], while the 
psedorapidity 77 = {l/2)\n[{p + pz) / {p — Pz)] = — ln[tan(0/2)]. For a relativistic particle 
{p ^ m), one gets y :^ r/ for 9 ^ m/p. 
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In order to get a non-zero integration region for Xi, one must take A> 1. It 
leads to 

r<xlcoshVut. (A. 12) 

The inequalities ( 1A.10|) and (]A.12p must be treated simultaneously with 



the inequalities (fT3|) . Depending on values of x±, the integration region splits 
into several parts: 

1. < a;_L < exp(— ?/cut)- In this case: 

x']_<T<x']_ cosh^?/cut , 
V^A-^ <xi< . (A. 13) 

2. exp(— ?/cut) < < (cosh?/cut)~^- Variables r, x\ run two subregions: 

< r < To , 

r < xi < 1 , (A.14) 

and 

To < r < x^ cosh^?/cut , 
^frA-^ < xi < VrA . (A.15) 

Here we introduced the notation]^ 

xx exp(-|/cut) 
^0 = 7^ r • (A. 16) 

2 - x_Lexp(?/cut) 

3. (cosh^cut)"^ < x_L < 1. In this case, variables r, x\ run the region ( fT3l) . 

Without cuts on rapiditites (formally, in the limit ?/cut — C)o), the full 
kinematically allowed region f|T3|) is restored. 



^Note that < tq < x\ cosh ^cut for the vahies of under consideration. 
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